Geometric properties of two-dimensional coarsening with weak 
disorder 
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Abstract. - The domain morphology of weakly disordered ferromagnets, quenched from the high- 
temperature phase to the low-temperature phase, is studied using numerical simulations. We find 
that the geometrical properties of the coarsening domain structure, e.g. the distributions of hull 
enclosed areas and domain perimeter lengths, are described by a scaling phenomenology in which 
the growing domain scale R{t) is the only relevant parameter. Furthermore, the scaling functions 
have forms identical to those of the corresponding pure system, extending the "super- universality" 
property previously noted for the pair correlation function. 



Introduction. — Many aspects of the out of equilib- 
rium relaxation of macroscopic systems still remain to be 
unveiled. The mechanism underlying coarsening phenom- 
ena, or the domain growth of two competing equilibrium 
phases after a quench from the disordered phase, is well 
understood [1] . However, an important part of the descrip- 
tion of this process is phenomenological and, to a certain 
extent, qualitative. 

A common feature of nearly all coarsening systems is 
that their dynamics are described by a dynamical scal- 
ing hypothesis with a single characteristic growing length 
scale R{t) that depends on the problem at hand. This 
hypothesis states that the domain morphology is statisti- 
cally the same at all times when lengths are measured in 
units of R{t). As a result, all time-dependences in correla- 
tion functions are encoded in R and lengths appear scaled 
by this typical length. The scaling hypothesis has been 
well-verified experimentally, numerically as well as ana- 
lytically in a small number of solvable simple cases. It is 
thus a well-established feature of phase ordering dynamics 
though no generic proof exists. 

Most studies of coarsening dynamics focus on the 
growth law R{t) and the scaling functions of various cor- 
relators and linear responses. An analytic determination 
of these functions remains open - apart from some toy 
models such as the one-dimensional ferromagnetic chain 
or the 0{N) model in the large N limit. A number of 



field-theoretic approaches to finding an approximate form 
of the scaling functions of two-point and two-time cor- 
relations have been proposed but none of them is fully 
successful [1]. 

In [2,3] we analysed the coarsening process in pure (not 
disordered) two-dimensional systems with non-conserved 
order parameter from a geometric point of view. We stud- 
ied the probability distributions of hull-enclosed and do- 
main areas and hull and domain- wall lengths. By deriv- 
ing an analytic expression for the number density of hull- 
enclosed areas we proved that the scaling hypothesis holds 
for this quantity. With a few additional assumptions we 
extended these results to the statistical properties of do- 
mains and boundary lengths and we put them all to the 
numerical test. In this Letter we extend these results to 
the coarsening dynamics of two-dimensional models with 
non-conserved order parameter under the effect of weak 
quenched disorder. By 'weak disorder' we mean random- 
ness that does not modify the character of the ordered 
phase (see [4] for a recent review). 

An important property of domain-growth in systems 
with weak disorder is the super-universality hypothesis 
that states that once the correct growing length scale is 
taken into account all scaling functions are independent 
of the disorder strength. This conjecture was first enun- 
ciated by Fisher and Huse [5] using a renormalization- 
group approach. The idea is that the length scale at 
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which the effects of quenched disorder are important is 
much smaller than the domain scale R{t). The latter 
dominates the elastic energy. Thus the dynamics of large 
structures is approximately curvature driven. Pinning at 
small scales modifies the scale factor, that is to say the 
growth law R{t), but not the scaling functions. The va- 
lidity of super-universality for the scaling function of the 
equal-time two-point correlation of several disordered sys- 
tems including the random bond Ising model was checked 
numerically in [6-9]. The stringest test proposed in [10] 
that includes higher order correlations was also passed nu- 
merically in [11]. 

In short, in this Letter we continue our geometric study 
of coarsening in two-dimensional systems. We simulate 
the dynamics of the random bond Ising model and we 
pay special attention to the scaling and super-universality 
properties of the probability distribution of areas and bor- 
ders. 

Hull enclosed area distribution. The ordering 
process of a magnetic system in its low temperature phase 
is visualized in terms of domains or regions of connected 
aligned spins. Each domain has one external perimeter 
which is called the hull. The hull enclosed area is the total 
area within this perimeter. In [2, 3] we derived an ana- 
lytic expression for the hull enclosed area distribution of 
pure, curvature driven, two-dimensional coarsening with 
non-conserved order parameter. Using a continuum de- 
scription in which the non conscrvc^d order parameter is 
a scalar field we found that the number of hull enclosed 
areas per unit area, n'l^{A,t) dA, with enclosed area in the 
interval {A, A + dA), is 



2ch 



{A + \htf 



(1) 



Ch = l/87r-\/3 is a universal constant that enters this ex- 
pression through the influence of the initial condition and 
was computed by Cardy and Ziff in their study of the 
geometry of critical structures in equilibrium [12]. \h is 
a material dependent constant relating the local velocity 
of an interface and its local curvature in the Allen- Cahn 
equation, v = — (A/i/2'7r) k [13]. The expression (1) is vahd 
for all quenches from equilibrium in the high temperature 
phase, Tq > Tc, to zero working temperature, T = 0. For 
a critical initial condition, Tq = Tc, the same expression 
holds with 2cfi replaced by Ch- Equation (1) can be recast 
in the scaling form 



nl{A,t) = 



{Xhtf ^ \\ht 



(2) 



with f(x) 



2ch 



In this way, scaling with the charac- 
^\ht, for coarsening dynam- 



teristic length scale, Rp{t) 
ics with non conserved order parameter in a pure system 

is demonstrated. The effects of a finite working tempera- 
ture are fully encoded in the temperature dependence of 



the parameter A/i while the same scaling function f{x) de- 
scribes n\ [3] as suggested by the zero temperature fixed 
point scenario [1]. 

When quenched disorder is introduced the growing phe- 
nomenon is no longer fully curvature driven and domain- 
wall pinning by disorder becomes relevant. At early times, 
the system avoids pinning and evolves like in the pure case. 
Later, barriers pin the domain walls and the system gets 
trapped in metastable states from which it can escape only 
by thermal activation over the corresponding free-energy 
barriers. In spite of these differences, coarsening in fer- 
romagnetic systems with quenched disorder also satisfies 
dynamic scaling [6-11] and a single characteristic length, 
R{t, T, e), can be identified (e is a measure of the disorder 
amplitude). As a result of the competition between the 
curvature driven mechanism and pinning by disorder, the 
coarsening process is slowed down and the characteristic 
radius of the domains depends on the disorder strength 
and it is smaller than the pure one, R{t, T, e) < Rp{t, T). 
Moreover, the super- universality hypothesis applied to the 
scaling function of the equal-time two-point correlation 
function in random ferromagnets was verified numerically 
in a number of works [6, 7]. 

The scaling and super-universality hypotheses suggest 
that Eq. (2) remains valid with the same scaling function 
f{x) = 2ch/{l +xf and {Xhtf^ replaced by R{t,T,e) 
for all 2d non-conserved order parameter coarsening pro- 
cesses in which the low-temperature ordered phase is not 
modified. More precisely, we expect that 



A 



(3) 



nh{A,t,T,s) = R-\t,T,e) / ^^2^^ ^ e 

should be valid in all these cases. In this Letter we test this 
hypothesis by following the dynamic evolution of the two- 
dimensional random bond Ising model {2d RBIM) defined 
by the Hamiltonian, 



H 



{hi) 



(4) 



where the Jij are random variables uniformly distributed 

over the interval [1 - e/2, 1 + e/2] with < e < 2. This 
model has a second-order phase transition between a high- 
temperature paramagnetic phase and a low-temperature 
ferromagnetic phase. We simulate the dynamic evolution 
of a model defined on a square lattice with linear size 
L = 10^, using a single-flip Monte Carlo technique and the 
heat-bath algorithm. Data are averaged over 10'^ samples. 
We show results for a random initial condition, crj = ±1 
with equal probability to mimic an infinite temperature 
equilibrium state, Tq — > oo. At the initial time t = 
we set the working temperature to a low value and we 
follow the evolution thereafter. We are interested on the 
effect of quenched randomness and we thus focus on a 
single working temperature, T = 0.4, at which the ordered 
equilibrium phase is ferromagnetic for all < e < 2. In 
what follows we drop the T dependence from R and n/j, 
so we simply denote them R{t,e) and nh{A,t,e). 
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Fig. 1: (Colour online.) Log-log plot of the characteristic do- 
main radius extracted from the collapse of the equal-time two- 
point correlation function against time after the quench. Dif- 
ferent curves correspond to different values of the quenched 
disorder strengths. For comparison we include the not disor- 
dered limit law R{t, £ = 0) oc y^. 

We determine the growth law, R{t,e), from a direct 
measure of the spatial correlation function C(r, t, s): 

1 ^ 

C{r,t,e) EE — ^(s,(t)sj.(t))||-_-.|^^ (5) 

i=l 

which also obeys dynamical scaling, 

C(M,e)^^^(T,s)/(^) , (6) 

with m(T, e) the equilibrium magnetisation density. In 
Fig. 1 we show R(t, e) as a function of t for several val- 
ues of the quenched disorder strength along with the pure 
case. We extracted R{t,e) by collapsing the equal-time 
correlations data using Eq. 6. After a few time-steps all 
curves deviate from the pure R{t, e = 0) oc ^/t law and the 
subsequent growth is the slower the stronger the disorder 
strength e. For strong disorder the characteristic length 
reaches relatively modest values during the simulation in- 
terval, for instance, R{t = 1024, e — 2) ^ 7 (measured in 
units of the lattice spacing). It is then quite hard to de- 
termine the actual functional law describing the late time 
evolution of R(t, e) and it comes as no surprise that this is- 
sue has been the matter of debate recently [14-21]. Using 
arguments based on the energetics of domain-wall pinning 
Huse and Henley [14] proposed the law R{t) oc (Inty/^ 
with ip = 1/4: for the random bond Ising model. More 
recently, powerful Montecarlo simulations [15-21] suggest 
a power law R{t) oc with a exponent 6 that depends 
on T and e. We are not concerned here with predicting 
the time-dependence of R. Instead, we use the numerical 
values in Fig. 1 to scale our data for the distributions of 
hull enclosed areas and hull lengths as explained below. 

In Fig. 2 we display data for one disorder strength, 
e = 2, taken at several times after the quench. Using the 




Fig. 2: (Colour online.) Number density of hull enclosed areas 
for one value of the disorder strength, e = 2, at several times. 
The data are shown in the scaling form and they collapse for 
A/R^ ^ 1. 

scaling form the data collapse for A/R"^ ^ 1 over 8 decades 
in the vertical axis and 4 decades in the horizontal axis. 
Deviations are seen for small areas and we discuss their 
possible origin below. This analysis confirms the scaling 
hypothesis in the region of large areas, A/R"^ ^ 1. 

In Fig. 3 we test the super-universality hypothesis as ap- 
plied to the hull-enclosed area distribution by presenting 
data for a single instant, t = 256 MCs, and several values 
of the disorder strength, e = 0, 0.5, 1, 1.5, 2. Again, for 
areas such that A/R^ ^ 1 there is a very accurate data 
collapse while for smaller areas deviations, that we discuss 
below, are visible. The solid black line represents the an- 
alytic prediction for the pure case that yields, under the 
super-universality hypothesis, the analytic prediction for 
all e. The curve falls well on top of the data. In the tail of 
the distribution we see downward deviations that are due 
to finite size effects already discussed in detail in [3]. 

Let us now examine several possible origins for the devi- 
ation of the numerical data from the analytic prediction at 
small areas, A/R^ ^ 1, in Figs. 2 and 3. The first source 
of problems could be the fact that we need to use a finite 
working temperature in the disordered case to depin the 
walls while the analytic results are derived at T = 0. In 
the pure case the effect of temperature is to create ther- 
mal domains within the genuinely coarsening ones. Based 
on this observation in [3] we explained that the small area 
probability distribution may have an excess contribution 
coming from these thermal fluctuations. By extracting the 
contribution of the equilibrium distribution of thermal do- 
mains, that became important for temperatures T ^ 1, we 
showed that the number density of hull-enclosed areas at 
T > is given by the zero temperature result once scaled 
by R{T,t). In the present disordered case, the working 
temperature we use is too low to generate any thermal 
domains and thus temperature cannot be the source of 
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Fig. 3: (Colour online.) Number density of hull enclosed areas 
at t — 256 MCs and several values of the disorder strength. 
The collapse of all data for A/R^ ~ 1 supports the validity of 
super-universality. The black solid line represents the analytic 
result for the pure case. 



deviations from the analytic form. 

Another possible origin of the difference between numer- 
ical data and theoretical prediction is the fact that the an- 
alytic results are derived using a continuum field-theoretic 
description of coarsening while numerical simulations are 
done on a lattice. In the presence of quenched random- 
ness one can expect the effects of the lattice discretization 
to be more important than in the pure case, especially 
for relatively small structures. Moreover, disorder induces 
domain- wall roughening and the AUen-Cahn flat interface 
assumption is not as well justified. 

Finally, the super-universality hypothesis is argued for 
large structures only [5] and thus the small hull enclosed 
areas are not really forced to follow it strictly when A ^ 

Geometric structure of the domains. In order 
to better characterise the geometric structure of the coars- 
ening process we study the relation between hull-enclosed 
areas and perimeters. In Fig. 4 we present the aver- 
age of the scatter plot of the scaled hull enclosed area, 
A/R'^, against the corresponding scaled perimeter, p/R, 
in a log-log plot for several values of the disorder strength 
at t = 256 MCs. We observe a separation into two regimes 
at the breaking point p*/R ^ 10 and A* / R^ ^ 2. In the 
upper and lower regime areas and perimeters are related 
by two power laws that do not depend on the disorder 
strength. More precisely [3], 



A 



R^{t,e) 



P 



R{t,e) 



(7) 



with 



Y> ^a> ^ 1.12 ±0.10 



for 



P 



1.83 ±0.10 for 



R{t,s) 
P 

R{t,e) 



10 



10 .(8) 



The upper exponent also characterises the highly ramified 
structures of the initial condition. Note that the small 
value of is not related to the existence of holes in 
these large structures, as suggested for domains in [22], 
since hull-enclosed areas do not have holes in them. The 
crossover between upper and lower regimes depends on 
time when observed in absolute value. In other words, dur- 
ing the coarsening process a caracteristic scale p* oc R{t, e) 
develops such that hull-enclosed areas with perimeter p > 
p* have the same exponent ao ~ 1.12 as in the initial con- 
dition before the quench [3] while hull-enclosed areas with 
smaller perimeter are more compact as indicated by the 
smaller value of the exponent a^. 
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Fig. 4: (Colour online.) Averaged scatter plot of the scaled 
areas against perimeters for several values of the strength of 
the quenched randomness. 

Hull length distribution. We now examine the 
hull length distribution that is expected to follow the scal- 
ing and super- universality hypotheses. In Fig. 5 we dis- 
play the number density of hull lengths a,t t = 256 MCs 
and e = 0, 0.5, 1, 1.5, 2. We show the data in the form 
suggested by the scaling hypothesis 



nhip,t,e) = R ^{t,e) g 



P 



(9) 



R{t,e) 

The data collapse on a master curve for all lengths sat- 
isfying p/R ^ 2 a value that is roughly the location of 
the maximum in the distribution. The super-universality 
hypothesis holds in this regime of lengths. 

Relating the distribution functions of areas and perime- 
ters through a change of variables, as explained in [3], we 
find 



R^{t, e) nh{p,t,e) 



a^Ch 



1 + 



R{t,e) 



(10) 



The scaling function is thus the same as in the pure sys- 
tem [3] with two branches characterized by the exponents 
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Fig. 5: (Colour online.) Number density of hull lengths at 
t = 256 MCs and several values of the disorder strength given 
in the key. The collapse of all data for p/R ~ 2 supports the 
validity of super-universality in this regime of lengths. 



and a^. As in the pure case the maximum in the 
distribution is described by Eq. (10) with a< while the 
maximum found with falls outside its range of valid- 
ity. 

In Fig. 6 we display the number density of hull lengths 
scaled with the typical radius R for one value of the disor- 
der strength at several times after the quench. The upper 
and lower predictions are shown with solid black lines. 
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universality for structures with A/R^ ^ 1 and p/R ^ 2. 

We showed that the analytic prediction for the num- 
ber density of hull enclosed areas derived for pure systems 
also describes the statistics of these quantities in the pres- 
ence of quenched ferromagnetic disorder. The geometric 
properties of the boundaries between phases are, in prin- 
ciple, more sentitive to quenched randomness than their 
interior. We showed, however, that the relation between 
areas and interfaces and, in consequence, the distribution 
of hull lengths are independent of the disorder strength 
also satisfying super-universality. 

In previous work all numerical tests of the super- 
universality hypothesis in the RBIM have focused on 
the study the equal time two-point correlation func- 
tion [6-9, 11]. The geometric approach that we have pre- 
sented enables us to show that this property applies even 
to very small structures. 

We also studied a number of related problems on which 
we report below. 

• We verified that analogous results are obtained for 
different probability distributions of the coupling 
strengths as long as these remain ferromagnetic. 

• We observed that the scaling plots do not depend on 
the working temperature while the latter is below the 
critical point. 

• All our results can be easily extended to the study 
of domain areas and walls as done in [3] for the pure 
case. 

• We checked that the super-universality hypothesis is 
also verified in c? = 1 [3, 23] for relatively long do- 
mains. Small deviations for very short ones are also 
observed in this case. 
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Fig. 6: (Colour online.) Scaling plot of the number density 
of hull lengths for a disorder with strength e = 1.0 at differ- 
ent times. The solid black lines are the analytic prediction in 
Eq. (10) and explained in the text. The part around the maxi- 
mum is characterized by the exponent while the tail of the 
distribution is well-described by the exponent . 



Conclusions. — In this Letter we analysed the statis- 
tics of hull-enclosed areas and hull lengths during the 
coarsening dynamics of the 2d RBIM with a uniform distri- 
bution of coupling strengths. We found that the number 
densities of these observables satisfy scaling and super- 
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